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Abstract 

We comment on the uniqueness of t-evolution(t = log{Q^/A^)) of non¬ 
singlet structure functions at low x obtained fromDGLAP equations 
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In recent years,an approximate method of solving DGLAP equations |^, 
ii at low X has been pursued In that approach, we expressed the 

DGLAP equations as partial diffetential equations in x (the Bjorken variable) 
X = Q‘^f2p.q and t{t = InQ "^/using the Taylor series expansion and assuming 
its validity at low x. One of the limitations of the approach is that the solutions 
reported are not unique]^, They are selected as the simplest ones with a single 
boundary condition-the non-perturbative x distribution at some initial point 
t = tQ. However, complete solution of DGLAP equations with two differential 
variables in general need two boundary conditions |^. 

The aim of the present note is to explore the uniqueness of the solution 
when it satisfies some given physically appropriate boundary conditions. The 
DGLAP equation for non-singlet structurefunction which evolve independent of 
singlet and gluon distributions [Q, |[ |[ Q is 

= [4^][3+41og(l-x)F^^(a;,t)+2 f 

at t l — z z 

(1) 

where t = log(^) and Aj = being the number of quark flavours. 

Let us introduce the variable u = 1 — z and note that 


OO 



( 2 ) 


This series (2) is convergent for u < 1.Since x < z < l,so 0 < u < 1 — a; and 
hence the converegence condition is satisfied.Using (2),we write in (1) 


pNS^ 


-,t) = F 
z 


NS 


OO / 




1^1 




dx^ 


(3) 


which covers the whole range of it, 0 < it < 1 — a;. Non-singlet structure function 
are expected to be well behaved in the entire x range,unlike the gluon or singlet 
structure functions which might diverge for a; —> 0 as in BFKL inspired models |^, 
p^.It is therefore justif ied if the higher order derivatives ^q^i for I > 1 are 
neglected in (3).This is more justifiable for small x{x 1),yielding 


, t) = t)+f; 

Z “ OX 




Putting (4) in (1) and performing the u-integration,one gets 

dF^^{x^t) dF^^{x,t) i^a 

t) --- - + P{x) -^-= R{x, t, 


where 


dt dx 

P{x) = -A/a;[2 log(l/a:) + (1 - a;)^] 


(4) 

(5) 

( 6 ) 
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and 


Q{x,t) = t 


( 7 ) 


R{x, t, = R'{x)F^^ix, t) (8) 

with 


R'{x) = -Af[3 + 41og(l - x) + {x - l)(x + 3)]. (9) 

The general solution of (5) is obtained by solving the following auxiliary 
systems of ordinary differential equations 

dx dt dF^^{x,t) 

P{x) Q{t) R{x,t, F^^{x,t) 


. If 

and 


u{x,t,F^^) = Cl 


vix,t,F^^) = C2 


are two independent solutions of (10) , then in general,solution of(5) is 


( 11 ) 

( 12 ) 


f{u,v)=0 (13) 

where / is an arbitrary function ofu and v. 

The auxiliary system (10) has three equations: 


II 

(14) 

dx dF^^ 

(15) 

P{x) R{x,t,F^S) 

and 


dt dF^^ 

(16) 

Q{t) R{x,t,F^sy 

Solving (14) and (15) one gets 


u{x,t,F^^)=tX^^{x) 

(17) 

v{x,t,F^^) =F^^{x,t)Y^^{x). 

(18) 


Solution of (16) is in general not possible. It needs additional informationof 
explicit X and t dependence ofF^^{x,t). In (17) and {18),X^^{x) andy'^‘®(x) 
are defined by 
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y-(x)=exp[/i?^l^ (20) 

Explicit analytic forms of andy^‘^(a:) in the leading l/x approximation 

are 

X^^{x) exp[-t log I loga:|] (21) 

while 

Y^^{x) w 1. (22) 

Note that (21) is strictly valid for 

— 11 — 

a:«exp[ ^ ^ (23) 

The possible form of (13) linear in is 

u + av = P (24) 


where a and /3 are constants to be determined from physically appropriate 
boundary conditions. 

In our earlier communications ^we set /3 in(24) = 0. 

The physically plausible boundary conditions for non-singlet structure func¬ 
tions are 

F^^{x,t) = F^^ix,to) (25) 

for some low t = and 

F^^{l,t)=Q (26) 

for any t. While the first one corresponds to a non perturbative input at some 
low momentum transfer,the second one corresponds to the expected largea; be¬ 
haviour of any structure function(singlet as well as nonsinglet) at any scale of 
momentum transfer consisten t within constituent counting rules SHI- Using 
the boundary condition of (25) and (26) in (24), we have 

toX^^ix) + aF^^ix, to)Y^^{x) = /3 (27) 


and 

tX^^{l) = p 

which leads to 


F^^{x,t) 


il)F^^ix,to) 

to 


^^S(1) _ ^JVS(^) 
AXWS(I) _ ■ 


As from (21) 


X^^{1) « 0, 


(28) 


(29) 


(30) 
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(29) yields 


( 31 ) 


to 

which was exactly our previous results ||, Thus,F^® ~ t behaviour at low x 
follows from the two boundary conditions (25) and (26) rather than the adhoc 
assumption /3 = 0 in (24). The auxiliary equation (14) can also have a solution^ 
u~^{x, t) ins lead ofit(a:, t) defined in (17).In that case, using the same procedure 
as earlier, we will have 


F^^{x,t) 


{l)F^^{x,to) 

to 


X{l)-^-X{x)-^ 
*fX{l)-t - X{x)-^ 


(32) 


which is ill dehned due to(30) in a; ^ 1 limit and hence excluded on physical 
grounds. 

To conclude,we have shown that the linear rise of F^^{x,t) witht at small 
a; is a unique prediction of the approximated DGLAP equation with physically 
plausible boundary conditions (25) and (26).Our result is howevervalid at very 
small X when condition (23) is satished.For its quantitative estimate,we define 
a ratio 


h{x) 


X 

exp[-i(l - x)2] ■ 


(33) 


As an illustration,it yields h{x) ~ 10“^ fora; ~ 6.10“^.This smalla; value is 
wellwithin the HERA regime ||l^ of a; > 6.10“®,but outside the available CCER 
neutrino data Q,where x > 7.5.10“® 
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